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Cosmologies with a Chaplygin gas have recently been explored with the objective of explaining 
the transition from a dust dominated epoch towards an accelerating expansion stage. We consider 
the hypothesis that the transition to the accelerated period involves a quantum mechanical process. 
Three physically admissible cases are possible. In particular, we identify a minisuperspace configura- 
tion with two Lorentzian sectors, separated by a classically forbidden region. The Hartle-Hawking 
and Vilenkin wave functions are computed, together with the transition amplitudes towards the 
accelerating epoch. Furthermore, it is found that for specific initial conditions, the parameters 
characterizing the generalized Chaplygin gas become related through an expression involving an 
integer n. We also introduce a phenomenological association between some brane- world scenarios 
and a FRW minisuperspace cosmology with a generalized Chaplygin gas. The aim is to promote a 
discussion and subsequent research on the quantum creation of brane cosmologies from such a per- 
spective. Results suggest that the brane tension would become related with generalized Chaplygin 
gas parameters through another expression involving an integer. 

PACS numbers: 98.80.-k, 98.80.Qc PU-ICG-04/08, |gr-qc/0404111| 



I. INTRODUCTION 



Mounting evidence from supernova la (SNIa) obser- 
vations U and the more recent cosmic microwave back- 
ground radiation (CMBR) data (see, e.g., Ref. @] and @) 
is suggesting that the expansion of our universe seems to 
be in an accelerated state. This has been designated as 
the "dark energy" effect Q. A possible candidate re- 
sponsible for this evolution is the usual vacuum energy 
represented by a cosmological constant, A, providing a 
negative pressure However, this idea has to deal 

with the fact that the observational value of A is 120 
orders of magnitude smaller than that established from 
field theory methods [U,E1- Moreover, it is a special "cos- 
mic coincidence" that the energy density associated with 
the cosmological constant has a value near to the density 
of other matter types exactly today. A plausible alterna- 
tive is to consider a dynamical vacuum energy (also called 
"quintessence") 0, involving one or two scalar fields, 
where some scenarios are associated with potentials jus- 
tified from supergravity theories Q . Nevertheless, when 
quintessence models confront the cosmic coincidence is- 
sue, they face fine-tuning problems and no satisfactory 
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solution has yet been found. 

Recently, an interesting proposal has been made Q, 
describing a transition from a universe filled with dust- 
like matter to an accelerating expanding stage: the sce- 
nario of the Chaplygin gas ap plie d to cosmology, which 
was later generalized in Ref. jiol [TT| . The generalized 
Chaplygin gas model is described by a perfect fluid obey- 
ing an exotic equation of state [111 ] 



A 

V = — 
P 



(1) 



where A is a positive constant and < a < 1. The 
standard Chaplygin gas corresponds to a — 1 . Insert- 
ing this equation into the relativistic energy conservation 
equation, leads to an energy density evolving as 



p{a) = A 



D 



a 3(l+a) 



(2) 



where B is a constant. This model interpolates be- 
tween a universe dominated by dust and a DeSitter stage, 
through a phase described by an equation of state of the 
form p — ap. It is curious to mention that the Chaplygin 
gas was originally introduced in aerodynamics |l2j |. 

Quite a few publications devoted to Chaplygin gas cos- 
mological models have alre ady appeared in the literature 
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Recent reviews can be found in |25l l2a | and a dy- 
namical system analysis regarding the stability of FRW 
cosmologies with a generalized Chaplygin gas was pre- 
sented in Ref. |27| . 

Within an observational context, it is worth mention- 
ing the analysis of compatibility at the level of CMBR 
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[3111 EH, namely from recent WMAP data 0, to- 
gether with SNIa observations as well as the mass 
power spectrum of large scale structure 0, ^1 an d uni- 
fied dark matter scenarios 0] . Overall, some favourable 
indications have been drawn when confronting general- 
ized Chaplygin gas models with cosmological observa- 
tions. 

Within a theoretical point of view, Chaplygin space- 
time models have also received increased attention. In 
particular, they can be formally related to an effective 
description of a complex scalar field whose action can be 
written as a generalized Born-Infeld action 0, 0, , 
corresponding to a "perturbed" d— brane in a {d + 1, 1) 
spacetime [18|. In particular conditions, a generalized 
Chaplygin gas can be obtained in the context of a general- 
ized Born-Infeld phantom framework [30l | . For a = 1 the 
equation of state associated with the parametrization in- 
variant Nambu-Goto d— brane action in a (d+1, 1) space- 
time can be retrieved [l^- Moreover, this action leads, 
in the light-cone gauge parametrization, to the Galileo- 
invariant Chaplygin gas in a (d,l) spacetime and to the 
Poincare invariant Born-Infeld action in a (d, 1) space- 
time. For an interesting survey relating fluid mechanics 
and c?-branes see also Ref. pgj . 

The perspectives indicated in the previous paragraph 
suggest that the analysis of a space-time with a Chap- 
lygin gas matter content also conveys the opportunity 
to investigate brane-world physics phenomenologically. 
Brane-worlds constitute a promising framework to ex- 
plore cosmological effects (in view of the recent WMAP 
data about inflation Q, as pointed out and discussed 
in Ref. [IJ), together with explaining other fundamen- 
tal physical issues (e.g., the hierarchy problem In 
the published literature the following scenarios can be 
found regarding brane cosmologies: the Arkani-Hamed, 
Dimoupolos and Dvali (ADD) model the Randall- 
Sundrum (RSI and RSII) models (see, e.g., Ref. |3^] and 
[34|. respectively) and the ekpyrotic model [3E1. For a 
recent review on brane-worlds see Ref. |36l l37| . Perti- 
nent issues concerning brane-world cosmologies include 
[38J: (i) Early evolution of the universe, (ii) perturba- 
tions on the brane/bulk and spectrum of fluctuations, 
(iii) stability of extra dimensions, (iv) onset of infla- 
tion, (v) dark matter/energy problems, (vi) corrections 
to Einstein gravity and (vii) novel issues in geometry 
(e.g., non-commutative geometries). Besides the previ- 
ous topics, brane-world cosmology also has to address the 
fundamental question of the origin of our universe. In- 
teresting proposals, regarding the quantum cosmological 
creation of brane world models, were recently introduced 
and analysed in Ref. Ea, E3, E3, |4^| . 

Within the physical context herein described, the main 
purpose of this paper is to investigate if a quantum me- 
chanical analysis would provide complementary informa- 
tion regarding FRW cosmologies with a Chaplygin gas 
content. In particular, assuming that after the universe 
reached a stage where dust matter was dominant, the 
subsequent evolution towards an exponentially expand- 



ing epoch involved a quantum mechanical transition as- 
sociated with some remnant component of the original 
wave function of the universe. In other words, underlying 
the classically observed universe on large scales there is a 
quantum mechanical background, where the cosmological 
influence of a wave function vestige (i.e., a rapidly oscil- 
lating state with a particularly small amplitude) would 
emerge during the dust dominace epoch. This specific 
wave function remnant had been present through that 
stage of the universe evolution and could be considered 
as a robust component with respect to decoherence pro- 
cesses in the early universe (see Ref. 0, 0, US for de- 
tails). In this framework, some pertinent questions may 
be raised. Could such a quantum description determine 
important physical consequences regarding a cosmolog- 
ical Chaplygin gas scenario? Namely, determining the 
transition probabilities towards an accelerated stage and 
restricting the range of values for the Chaplygin gas pa- 
rameters? 



We also focused on an additional aim: a discussion on 
the quantum creation of a brane-world cosmology. This 
perspective is suggested by the phenomenological con- 
text that some brane-worlds acquire within generalized 
Chaplygin gas cosmologies. This allows us to make use 
of important similarities between the quantum cosmology 
of a brane-world model and a FRW quantum cosmology 
with a Chaplygin gas content, which we will point in 
detail herein. While the former may require some non- 
standard tools (see e.g., Ref. US), the latter imports 
techniques from well known quantum cosmological devel- 
opments more detail, an analysis can 
be made exploring the minisuperspace resemblances be- 
tween particular brane-worlds, generalized Chaplygin gas 
cosmological models and FRW minisuperspace scenarios 
with a complex scalar field cj), with (j> = e ze ip (where 9 con- 
stitutes a cyclical coordinate, determining that the con- 
jugate momentum irg corresponds to a conserved charge). 
We will present herein some interesting results which are 
retrieved when including additional quantum corrections 
to the WKB description. 



Accordingly, this paper is organized as follows. In 
Sec. II we present our minisuperspace model, pointing to 
the possible relation between a generalized Chaplygin gas 
cosmology and some brane-world properties. The cor- 
responding Hartle-Hawking and Vilenkin states are ob- 
tained in Sec. III. Interesting quantization consequences 
are described in Sec. IV, and possible extensions beyond 
the saddle-point approximation are discussed in Sec. V. 
Sec. VI constitutes a summary of the results herein pre- 
sented. An appendix is also included. It conveys the cal- 
culations of turning points, which are pertinent for the 
minisuperspace case (investigated throughout sections II 
to V), whose interesting similarities with some brane- 
world cosmologies are discussed in this paper. 
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II. MINISUPERSPACE MODEL 

Let us take the case of a closed (k = 1) FRW model 
in the presence of a positive cosmological constant A > 
with a generalized Chaplygin gas 1 . The minisuperspace 
Lagrangian has the form (see, e.g., ref. [531 154^ 



L 



3vr 
AG 



N 



Na - 2i:a 3 N' 



A 
8wG 



2ira 3 Np. (3) 



After determining the corresponding canonical momen- 
tum to a and therefore the Hamiltonian constraint, 



G 
"to? 



+ §V(a) = 0, 



(4) 



its quantum description can be represented by the 
Whecler-DeWitt equation, written as 2 



G d 2 3tt , , 



*(o) = 0, 



(•5) 



where ^(a) is the wave function characterizing the quan- 
tum state of this FRW universe. The effective minisu- 
perspace potential V(a) is given by 



V(a) 



Aa 4 — p(a)a , 



(6) 



where A 



In particular, ipx and ip2 constitute outgoing or ingoing 
modes in the classically allowed regions, decreasing or 
increasing functions in the classically forbidden regions, 
respectively, and C\,C<x are constants. The WKB valid- 
ity condition is given by G 



dV(a) 



da 



«|V(a)| 3/2 . 



In order to retrieve the wave functions in (JHJ) in an ana- 
lytical form, a saddle-point formulation will be employed. 
Therefore, it is necessary to determine the possible turn- 
ing points for the potential V(a). Analyticity also re- 
quires we employ the following approximation (where we 
henceforth drop the overbars on A and B) 



A 



B 



l + n 



o 3(X+c0 

1 1 



2 1 



B~° 
1 



1 A 

1 + a B 
A 2 



a 3( Q + l) + 



B 2 



- a 6(«+l) 



and if we take 



a < 



3( = + l) 



we can write the potential V(a) as 

V (a) w a 2 - Aa 4 - aB^ 



(11) 



(12) 



(13) 



p(a) = [A + 



B 



a 3(l+a) 



(7) 



together with < a < 1, A = A(^f-) a+1 



B = 

As mentioned, for a = 1 we 



B(8E£) a+ \ A,B e i 
recover the usual Chaplygin gas case. We will consider 
a semiclassical WKB (Wentzel, Kramers, Brillouin) de- 
scription, where the wave function of the universe can be 
approximated by 



where 



4>i 



exp 



S«(a) 



(8) 



(9) 



with i = 1,2 and Sq\o) satisfying the Hamilton- Jacobi 
equation 



'dS^(a) 



da 



9tt 2 



(10) 



1 See Sec. VI for a discussion on the presence of A: it makes our 
model analytically feasible and can represent the lowest order 
effect of the effective cosmological constant determined by the 
generalized Chaplygin gas. 

2 We are ignoring factor ordering aspects related with 7r 2 , with n a 
being the conjugate momentum to a. 



where we neglected the second as well as higher order 
terms in JTTJl as < a < 1 (sec Eq. © and JTTJ)). We can 
thus identify a suitable range of values for the parameters 
A, B, a, and the scale factor a where Eq. (|13fl is valid, 
corresponding to the analysis presented in this paper. 
As far as the next order of approximation is concerned, 
we will point in section V how the term proportional to 
a 4+3a j n y( a ) can be dealt with in a perturbative manner 
when A <C 1. 

In the linear approximation provided by Eq. <I13|) wc 
can identify three physically different cases |60l l61| . A 
plot for V~o(a) regarding the second and third cases is 
given in Fig. ^ (see appendix A for more details): 

1. When KB 2 > |, where B = B^, there is one 
real root and pair of complex conjugated roots for 
the equation /(a) = a 3 — a/A + B/X = 0. The 
only value for which Vo(a) vanishes is a = and 
V (a) < for a > 0. 

2. In this case, we have B — |-^- It is found that 
Vq (a) vanishes for a — and there is a positive root 
for f(a) = when a = 1/vA. The other root is at 
a = -2/VA. We further have that V Q (a) < for 
< a < 1/VX and a > l/VX. 

3. In the final case, we have AB 2 < | and three 
real roots can be found such that two are posi- 
tive: ax > a 3 > > a-i- Moreover, Vo(a) < 
for < a < 0,3 and a > a\, but Vo(a) > 



4 



V D (a) 



o.i 
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FIG. 1: This figure shows the behavior of the potential Voia) 
as function of a. The dashed line corresponds to the case 
AB 2 = 4/9. In this case, Vo{a) vanishes at a — l/v^A and 
there are two Lorentzian zones corresponding to a < 1/%/A 
and l/VX < a. The solid line corresponds to AB 2 < 4/9. 
The potential vanishes at a = ax and a — a?,. In this case 
there are two Lorentzian zones (a < a-i and a\ < a) separated 
by a classically forbidden region (a3 < a < ai). 



for d3 < a < a\. We thus have a classically al- 
lowed region followed by a classically forbidden Eu- 
clidean region and finally there is a second classi- 
cally allowed Lorentzian region. For < a < 
there is always a region such that the linear ap- 
proximation (|13fl is valid. For 03 < a < a\ 
this approximation is valid for a\ <C a m ax where 

«max = [(<* + l)f ]3("+d. For a > ax, if ai < a max , 
the linear approximation is valid if there are val- 
ues of a such that a\ < a <C a max . This leads to 

/ \3(a+l) , 

f-^=J < -f, which together with KB~^ < | 

implies that 3 3 < a+1 UB 2 < 1. 



Before proceeding to obtain quantum solutions for the 
Wheeler-DeWitt equation (see Eq. J^J and (fT^l l. it is of 
relevance to point the following. Quantum cosmologi- 
cal scenarios with a physical configuration similar to the 
third case above described (where AB 2 < | ), can be 
found in the literature but in somewhat differently re- 
lated contexts. 

One such situation is present in Ref. |55l I5d| when 
a FRW minisuperspace with a complex scalar field of 
the form <p = s te (p, with a mass term proportional to 
ip 2 and a conserved charge (i.e., conjugate momentum) 
associated with -Kg is considered. This produces a min- 
isuperspace potential U(a, tp; ttq) where ire is a constant 
parameter. For arbitrary constant values of ip we obtain 
different sections in the potential U, inducing a similar 
dependence on a as in Vo(a) given by Eq. (|13fl and de- 
picted in Fig. n Generalized tunnelling (Vilenkin) jE^ 
and Hartle-Hawking |5?1 | boundary conditions were then 



introduced in Ref . l55l l56l . The reason that the physi- 
cal model in Ref. |55ll56| relates to the above Chaplygin 
gas Hamiltonian formulation could be traced to the fact 
that the Chaplygin gas FRW action can be retrieved from 
the action of FRW models with the complex scalar field 
(j> = e l9 ip, under specific conditions [9l flOl fll| . 

The association between generalized Chaplygin gas 
cosmologies and minisuperspace models with complex 
scalar fields can be brought towards a recent and promis- 
ing domain: that of brane cosmology l3a.|3j|. In fact, the 
complex scalar field employed in Ref. [UJ, to obtain 
a generalized Chaplygin gas can be identified within a 
Born-Infeld action, which arise form the imbedding of the 
(3+l)-D brane into the (4+l)-D spacetime More- 
over, it can also be seen to correspond to a "perturbed" 
d— brane in a [d + 1, 1) spacetime 01- Furthermore ac- 
tions of the Born-Infeld have recently been the subject of 
wide interest (see ref. |2^| and references therein). This 
comes comes from the result that the effective action for 
the open string ending on d— branes can be written in a 
Born-Infeld form. In that respect, Born-Infeld cosmolog- 
ical solutions could assist us in the understanding of some 
brane-world dynamics regarding the universe evolution. 

The above context regarding brane-worlds can be fur- 
ther specified. In fact, FRW models with a Chaply- 
gin gas may constitute a phenomelogical tool to study 
some brane-world models: another physical situation 
with two Lorentzian sectors separated by a classically 
forbidden Euclidean region can be found in the Hamil- 
tonian treatment of a DeSitter (3 + 1)— brane imbed- 
ded in a (4 + 1)— Minkowski background [3|j. The pres- 
ence of a classically disconnected epoch (as in Fig. 1 for 
a < 03) has been identified as a property of brane-world 
quantum cosmology |39|. These remarks thus suggest 
the view that FRW models with a generalized Chaply- 
gin gas content would enable us to investigate specific 
aspects of particular brane-world models. In particular 
through an Hamiltonian formulation, where symmetries 
and conserved quantities, possibly related with physical 
properties of brane cosmologies, will be of relevance. In 
Sec. V, where we consider the next order of approxima- 
tion in Eq. © and (|13fl (going beyond the Hamilton- 
Jacobi equation (|10fl L we will point a further correspon- 
dence between brane-world parameters and a generalized 
Chaplygin gas. 



III. HARTLE-HAWKING AND VILENKIN 
WAVE FUNCTIONS 

We begin by addressing the second case (j3 = 

we have two classically allowed Lorentzian regions sepa- 
rated by the point a = l/VA where Vo(a) vanishes. Let 
us first analyze the region where a > 1/^/A. Besides 
having to satisfy the range present in Eq. (|12fl it is also 



5 



necessary to have 
1 



A 



(a+1) 



B 

A 



(14) 



(in order to have values of a such that the linear approx- 
imation is still valid) which leads to 



In the next paragraphs we will retrieve for the third 
case AB 2 < | the Hartle-Hawking jK^ and tunneling 
(Vilenkin) [Hsj wave functions (see also Ref. 0, Hil El 
l52j ] for more details). In order to obtain the solutions 
So (a) of Eq. (fTT?f) it is necessary to determine the analyt- 
ical expressions 



b 4 < a ; 



a+1 



(15) 



The wave function is then expressed in the form of Eq. 
and 10 where 



S (a) = i 



3tt 



x/-Vo(a)da. 



(16) 



\ "A 



Introducing the Lorentzian conformal time r] using 

dSo . 3tt da 
da 2 dr\ ' 

we get the following equation for a(rf) 

da 



(17) 



drj 



(18) 



from where it can be found that 



a (^) = —rr 



s/X Vcosh(?7) - 2 



(19) 



with r\ G ] — oo, — cosh -1 (2) [, noticing that when rj 
-co then a — > and when T] — > — cosh -1 (2) then a 
+oo. For the wave function it is obtained that 



vA 



^-V a (a)da = - 



a a + 



x/A 



a a + 



+ — In 

aI 



a a + 



VA 

\/3-ln (2 + \/3)/VA 



Ai 



VA 



(20) 



Employing the same method for the region where a < 
1/VA it is found |H3,E3 that 



VA V cosh(?7) + 2 



(21) 



with ?7 S [0, +oo [, noticing that when r\ — > +oo then a — > 
^= and when 77 — s- then a — > 0. For the wave function it 

can be concluded that J ^ ^/—Vo(a)da is equal to minus 
the expression in Eq. 



V\V (a)\da= x /- I I(a)da. (22) 



This integral will of course depend on the specific chosen 
region, e.g., between < a < a 3 for the first classically al- 
lowed Lorentzian region. Let us determine the values for 
the integral in ll'L'l) in the following. After some lengthy 
calculations, it is found that |6Cll61| 



ai a 2 {a - a{) 
I(a)da = y - P j 

/ ^2 2 m ( 2ai + a 3 
+ 7 |a 1 a3(a 1 -a 3 )n^,^- T ^ )g 

+ [a 1 a 2 (a 1 +a 3 )-6f3 2 +(3a 2 3 ]F{v,q) 



(ai + a 3 y 



2a 3 + 01 
+ (2 ai + a 3 )F(v,q)}} 



(a 3 - ai)n(w, l,g) 



(23) 



where /3 



~ 3 S 7 



a 3 + ai 



a3(2ai+a 3 ) 
ai(2a 3 +ai) 



u = arcsm 



(2a 3 +a 1 )(a-a 1 ) t thef w j th 
(2ai+o 3 )(a-a 3 ) ' ° 



I{a)da 



ai a 2 (a 3 -a) 

7 

7|aia 3 (ai + a 3 ) 







(a 3 - ai) n ^x, 9 ! + " ' /T ' ( \ • '/ 1 

(a% + a 3 f 



2ai + a 3 

(01 - a 3 ) n (x, g 2 , g) + (2a 3 + 01) F ( X , q) 
6 h 3 2 F( Xl q)+P~/a 2 F( Xl q)}, (24) 



with x 



arcsin -i / ai l a3 a \ , as well as f ai I(a)da 

V a 3 (ai— a) ' Ja 3 V / 
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lim a ^ ai J a a 3 1 (a) da where 
I(a)da 



a£ a 2 (a - a 3 ) 
3 P J 



+ 7 |— aiag (ai + a 3 ) II ^<5, — l — — -,r 
+ 6(3 2 F{S,r) 



/3a i 



( ai -a 3 )F(o,r) + a 3 n(o,l,r) 



+ /3(a x + a 3 ) 



a 3 n(6,r 2 ,r)-(2a 3 + ai )F(5,r)}} 



(25) 



with <5 = arcsin i / ? l( - a a , 3 ^ and r = . 

y (01—03)0 V (2a3+ai)oi 

The functions F and II are elliptic integrals of the 
first and third kind, respectively (see, e.g., Ref. [§3, 
l6ll | for more details). When calculating the expres- 
sion lim a _> 0l J I(a)da it is found that it is well de- 
fined except for the terms corresponding to a ( a ~ Q3 ) + 

\ Ia s a (i-a~°j^ a ' as eacn °^ tnese expressions diverges 
when a — ► ai. However, it can be seen that these diver- 
gences cancel each other |(3(3] and therefore the integral 
in (|25[1 is well defined and can be used to determine the 
transition amplitudes for the FRW universe with a gener- 
alized Chaplygin gas to tunnel from a Lorentzian region 
between < a < a 3 to another where a > ai (see Fig. 0. 

Following Ref. [55l IHfij the no-boundary (Hartle- 
Hawhing) wave function representing the region under 
the potential barrier Vo(a) corresponds to the increasing 
modes 



\P tx exp 



3tt 
2G 



a<ai 



\/Vb(a)da 



(26) 



where we are using Vq(cl) ~ — Aa(a — ai )(a — a 3 )(a + ai 
a 3 ) . In more precise terms wc obtain 59] 



*(a) = Ccos 
\E"(a) = Cexp 
*(a) = 2(7 cos 



3tt 
2G 



V-Vb(a)(ia + 



, (a < a 3 ) 



(27) 



(a 3 < a < ai) 



(28) 



3tt r 

2G 



y/-Vo(a)da - ^ 



(ai < a), 



(29) 

and C is an arbi- 



where C = Cexp §§ f^ 1 ^/Vo(a)da 
trary constant. 

As far as the tunnelling (Vilcnkin) wave function is 
concerned, it contains only outgoing modes and we find 



*(a) = Cjexp 



3tt 
2G 



X cos 



4z exp 



y/V (a)da 



X COS 



3tt 
2G 



V^W(a)da 



(a < a 3 ) 



*(a) = Cjexp 
+ 2i exp 



2g ^ V^Rda 



(30) 



, (a 3 < a < ai) 

(31) 



*(a) = 2Ce" r/4 exp 



3^ r a ' 

l 2G J ^~ V ^ a ) da 



(ai < a). 
(32) 



Having established explicitly the analytical form for 
the Hartle-Hawking and tunnelling (Vilcnkin) wave func- 
tions in Eq. if ^ -iP ty and Eq. pD jl -ip ^ . respectively, we 
can now analyze the transition amplitude A, for the uni- 
verse to evolve from the classically allowed Lorentzian 
region between < a < a 3 to the other where a > ai 
(see figure ^| 3 . This transition amplitude is given by 



A = exp (t2/) 



(33) 



where the — and + sign correspond, respectively, to 
the tunnelling (Vilenkin) and Hartle-Hawking wave func- 
tions, and we use the expression 



3ir f 1 

I=—J y/Xa(ai-a)(a 



+ ax + a 3 )(a — a 3 )da. (34) 



We henceforth consider the new dimensionless quantity 
XIG which will be analyzed as function of the parameter 



X = B^VX 



(35) 



that is also dimensionless. When X = we recover 
the well know closed FRW-DeSitter quantum cosmology 
case. From Eq. (|34|l it can be seen that XIG decreases 
when X varies from to 2/3 and hence the transition am- 
plitude, A, for the tunnelling (Vilenkin) boundary con- 
dition increases when X increases (for a fixed value of 
the cosmological constant), while for the Hartle-Hawking 
boundary conditions A decreases. 



3 For a discussion on the general meaning of tunneling universes, 
see, e.g., Ref. |62|. 
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In terms of the generalized Chaplygin gas parameters 
B and a, we have that when a increases from to 1 and 
B is constant and smaller that 1, X also increases and 
consequently XI G decreases, which for a fixed value of A 
determines that I also decreases. This then implies that 
the transition amplitude for the tunnelling (Vilenkin) 
boundary condition increases but the transition ampli- 
tude associated to the no-boundary proposal decreases. 
The opposite conclusion is reached when a increases and 
B is constant and larger than 1. Taking now fixed val- 
ues of a and A, if B increases then X increases. As 
a consequence, we conclude that similarly to the previ- 
ous situation (for a increasing when B is constant and 
smaller than 1) also the transition amplitude for the tun- 
nelling (Vilenkin) boundary condition increases but the 
transition amplitude associated to the Hartle-Hawking 
proposal decreases. In Fig. [21 we plot how XIG varies as 
function of X. 




FIG. 2: This figure shows the behavior of XIG, where / is 
defined in Eq. 1341 . in terms of the dimensionless parameter 
X defined in Eq. 13511 . For X = we recover the closed 
FRW-DeSitter case, while the case X — 2/3 corresponds to 
AB 2 = 4/9. In this last case there is no tunnelling between 
the two Lorentzian regions. 



latter we have 



R = AiA + 2nG 



A- 



B 



a 3(l+a) 



<L4 



B 



a 3(l+a) 



(36) 



As it can be shown, for small values of the scale fac- 
tor, the Ricci scalar R diverges. Furthermore, for small 
values of the scale factor, R in Eq. $jt$ will have the 
same behavior as that we would obtain from the approx- 
imations in Sec. II. Such divergence could be avoided if 
some other material content dominates the energy den- 
sity of the universe at that stage (for example some ef- 
fective scalar field), but we disregard this possibility in 
the present work. 

The presence of a curvature divergence becomes more 
complicated in the FRW scenario with a generalized 
Chaplygin gas: the divergence is present in a Lorentzian 
(classically allowed) region (see Ref. [3^ as well as 
[55l |H(|). A suggestion has been advanced, designated 
as DeWitt's argument [Hill Hi, in that #(a = 0) = 
is imposed. The interpretation is that the divergence is 
not replaced by an Euclidean conic-singularity-free pole, 
but instead is neutralized by making it quantum mechan- 
ically inaccessible to wave packets. 

Before proceeding, let us also mention the following. 
Taking A = in Eq. ©, it can be shown (see Sec. VI 
for more details) that V(a) for A <C 1 is essentially as 
represented in Fig. 1 for Vo(a) with A ^ (within the 
mentioned approximations). Hence, the validity of the 
results presented in this section concerning B and a are 
not restricted to our approximations. The physical con- 
sequences will hold in a wider scenario, albeit with some 
modifications in the mathematical expressions. 

Assuming therefore that the wave function vanishes at 
the origin 4 , i.e., a — => *(a = 0) = 0, from Eq. J2HJ), 
(|27|l it is obtained that in the case of the Hartle-Hawking 
wave function we will have 



n € Z. 



(37) 



IV. QUANTIZATION CONDITION INVOLVING 

B AND a 



Let us now establish how the quantum states obtained 
in the previous section may determine which sets of val- 
ues for B and a are allowed. 

In the case of a closed FRW universe with a positive 
cosmological constant and a generalized Chaplygin gas, it 
can be seen that a singularity is present at a — as the 
Ricci scalar diverges. This divergence in R is present 
whether we consider the approximation introduced in 
Sec. II [see Eq. lfTT|) - ((L3l ] or the the full expression of 
the energy density of the Chaplygin gas. In fact, for the 



It can be further seen that when B 1 + a is small 



^/-V (a)da ■ 



1 7T _ 2 



-a 2 + aB 1 *" = —B 1 +° 



From Eq. J^JJ and ® it follows that 
3tt 



AG 



Bi+a — 1 w An, 



(38) 
(39) 



with n G Z. 



4 See Sec. VI for a discussion on this issue. 
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Recalling that the Vilenkin wave function can be writ- 
ten as 



5" (a) = Cexp 



3tt 
2G 



y/V (a)da 



x cos 



3tt 
2G 



+ 4iGexp 
2G 



^-V a (a)da + - 



y/-V (a)da - 



(40) 



for the case of < a < and assuming again that 
\&(a = 0) = 0, it now results that 



_>(-; 



kit, 
kit, 



(41) 



where k, k G Z. Both conditions in l|41(l cannot gener- 
ally hold simultaneously. However, there is an exception. 
When |g f® 1 y / Vo(a)da ^> 1, it is enough that the last 
condition holds. This statement can be proven as follows. 
Let us take 



cos 



3^ r a 3 n 

2gL V=%W*,+ 4 



+4i exp 



cos 

If we now define 
T = 



3tt 
2G 



\/-Vn(a)da - 



= 0. 



exp 



3?r 

~g 



yfVo(a)da 



we can write 



1-4T 
1+4T ' 



(42) 



(43) 



where n € Z. These equalities represent l|41|) and then 
cannot both hold unless T ^ or T > 1. As we are 
considering the possibility that the universe tunnels from 
one Lorentzian region (a < 03) to another Lorentzian re- 
gion (ai < a) , the parameter T is greater that the unity. 
Consequently, the case T — + is physically impossible, 
while the case T 3> 1 is admissible. In this situation, we 
have, with h = 2k + 1, that 



3tt 
2G 



\J~V {a)da 



3tt 
T 



fc7T, 



(44) 



which is equivalent to the second equality in expression 

ffil; 

Finally, we can address the question when is the condi- 
tion |g f* 1 yJVo(a)da ^> 1 fulfilled. To this aim we write 



We know the behavior of XIG in terms of the parameter 
X [defined in Eq. (|35|l ]: XIG will take its largest values 
when X is small. So in order that exp[(37r/2)/] ^S> 1, it 
is necessary that XIG gets its largest values and AG be 
small, i.e., A has to be small. Hence, what we have is that 
the condition ^(a = 0) = applied to the Vilenkin wave 
function implies the universe to be close to a DeSitter 
model and to have a small cosmological constant. In this 
case we obtain that 



3tt 
4G 



B^ 



Ak. 



(46) 



V. BEYOND THE HAMILTON- J ACOBI 
EQUATION 

In this section we will consider the next order of 
approximation in Eq. (|llfl . i.e., bringing the terms 
a 4+3 "(l + a) -1 AB~~^ into the minisuperspace poten- 
tial. This will be done perturbatively 5 with the parame- 
ter A satisfying A <C 1, allowing a quantum cosmological 
analysis beyond the WKB formulation in Eq. (|5|)- (|10fl . 
Furthermore, it will allow us to include explicitly in the 
analysis a parameter that can be related to the brane 
tension in some brane world cosmologies (see Ref. 0]). 

Let us consider again the Wheeler-DeWitt equation 
(jSJ) and employ the following expansion where we take 
A< 1: 



tf(o) 



exp 



G 



(S (a) + AS 1 (a) + A 2 S 2 (a) + ---) 



(47) 



Using the expansion in Eq. l|l l[l up to the term in 
we then obtain in powers of A 



,4+3a 



1 



3ttG 



1 

3^ 

3tt 



dS {a) ^ dSija) 
da da 



v d 2 S (a) , .d 2 S 1 (a) 



da 2 



da 2 



A 2 dS 2 {a) 
da 

A 2 d 2 S 2 (a) 
da 2 



1 2 



*(a) 
*(o) 



+ 4G ^ o(a) + AVl ^ + A2y2(a) ] * (a) 



where 



V (a) 



1 2 - Aa 4 - aB-rh, 



Vi(a) = -- 



,4+3a 



-B~ — 



0, 

(48) 

(49) 
(50) 



with V2 (a) following in a similar manner (but will not be 
necessary at this level). 



exp 



3tt 
2G 



y/V {a)da 



exp 



2 



5 See Ref. loll and lo'J for a formal description regarding the ex- 
(45) pansion of a wave function of the universe in terms of a pertur- 

bative parameter. 
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In the WKB approximation where d f°i a ' > is neglected, 
we get at the A and A 1 level, respectively, the equations 

'dSo{a,y 



da 

2 dSo(a) dS\ (a) 



3nG da 



da 



9tt 2 

1 d 2 S 1 {a) 
3n da 2 



(51) 



gMa)=0, 



(52) 



where the former is the Hamilton- Jacobi equation (ob- 
tained in Sec. II). If we further assume that Si(a) is 
slowly varying, then Eq. I|52[l can be approximated by 



2 dSo(a) dS\(a) 
3ttG da da 



(53) 



If we consider that the turning points are exclusively 
determined by the approximate potential Vq and that V\ 
is a correction when A is small, then it can be possible to 
obtain (by an approximate analytical expression) a quan- 
tification rule involving a, B as well as A. In order to get 
analytical expressions for this case, we should integrate 
the expression 



J 



Vi(a) 



o V^Ma) 



zda, 



(54) 



as implied from Eq. (|51|l and l|53[l . By means of an ap- 
proximation similar to the one given in Eq. 1)38(1 , namely 
with iJ 1 /( 1 + Q ) being small (i.e., KB 2 small), then J can 
be approximated as follows 



(£?i+o _ a ) 2 da. 



This expression can be written analytically as 



(55) 



(56) 



where B is the beta function (see Ref. |s3,E3)- 

tification rules follow when using wave functions in the 
region < a < a 3 for the two boundary conditions that 
we have considered. We will not pursue this issue further 
in this paper, as it will be addressed in a future work. 

Therefore and with suitable choices of boundary con- 
ditions we can then obtain an explicit form for Si (a) and 
hence for 



W(a) w exp 



^(So(a) 



■AS^a)) 



(57) 



Within this procedure we can thus include explicitly the 
parameter A in the wave function of the universe. This 
may be of relevance: a connection between the Chap- 
lygin gas (a — I) and brane- world models is that the 
brane tension is In the case of a generalized 

Chaplygin gas the "generalized" tension would be 



(see Ref. 0). Therefore, if Si (a) is computed then it 
would also be possible to extract a quantization condi- 
tion involving A together with B and a, broadening the 
scope of Eq. (gSJl and g^). 

Accepting that the minisuperspace similarities pointed 
out in Sec. II between some brane- world cosmologies and 
FRW models with complex scalar fields (inducing the 
generalized Chaplygin gas matter content) provide a use- 
ful association, then these results may be of interest. Of 
course the range of models herein considered is yet re- 
stricted, but it may direct to potential insights regarding 
a quantum mechanical description of brane-worlds to- 
gether with the creation of corresponding cosmological 
models. 



VI. CONCLUSIONS AND DISCUSSIONS 

As a summary, let us point the main results conveyed 
in this paper. Our purpose was twofold. 

First, it was our aim to present a quantum mechan- 
ical description of a closed FRW model with a positive 
cosmological constant and a generalized Chaplygin gas. 
The physical scenario we are contemplating is that, un- 
derlying a classical universe, there was a wave function 
vestige (i.e., a rapidly oscillating state with a particu- 
larly small amplitude) whose physical implications would 
emerge within the epoch of cosmological dust dominance. 
Such specific remnant component of the original wave 
function of the universe had been present during that 
period of the universe evolution. It can be interpreted 
as a robust component (with respect to decoherence) in 
the sense introduced in Ref. [4^, |43, Hlg • In other words, 
after the universe reached a dust stage, the transition to- 
wards an accelerating expansion would be associated to 
a quantum mechanical transition for that specific wave 
function remnant component of the original state. In 
this context, we obtained the Hamiltonian and the corre- 
sponding Wheeler-DeWitt equation associated with such 
transition. In particular, we established a range for the 
variables and parameters of the model where a linear 
approximation for the effective minisuperspace potential 
could be used. 

The quantum mechanical amplitudes for the transi- 
tions from a dust stage towards an expanding accelerat- 
ing phase were obtained either in the form of a Hartle- 
Hawking or Vilenkin wave function. We established how 
these transition amplitudes vary according to the allowed 
values for a and B: when either a or B increase, the 
Vilenkin state (generally, see Sec. Ill) allowed the prob- 
ability for the universe to tunnel towards the accelerated 
stage to increase, while the opposite situation occurred 
for the Hartle-Hawking wave function. 

Subsequently, we investigated whether this quantum 
mechanical formulation for the transition from a dust 
stage towards an exponentially accelerated phase deter- 
mined any physical consequences regarding the Chaply- 
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gin gas cosmological scenario. We then found that the 
Hartle-Hawking and Vilenkin wave functions could im- 
ply that the parameters a and B would be related and 
restricted to a quantization condition 6 characterized by 
the presence of an integer n. 

This result is obtained as we impose the condition 
^f(a = 0) = corresponding to the presence of a di- 
vergence in the curvature. The divergence is generic, i.e., 
it is not a consequence of our linear approximation. The 
presence of such divergence in a classically allowed re- 
gion can be dealt following, e.g., the approach suggest 
in Ref. [3^, |5^, |5(| (also designated as DeWitt's argu- 
ment): the divergence is not replaced by an Euclidean 
conic-singulatity-free pole, being neutralized by making 
it quantum mechanically inaccessible. Let us remember 
that, e.g., the tunneling (Vilenkin) boundary condition 
was thought of for a problem of (semi-classical) creation 
of the universe, in analogy with the tunneling process in 
quantum mechanics. Note, however, that the singularity 
in our model is located in a classical region, before tun- 
neling. We can consider waves proceeding from a classi- 
cal region, with a probability to exit from it towards the 
other classical region (see figure 1). 

Nevertheless 7 , it is only a sufficient but not necessary 
that the wave function vanishes at the origin in order to 
avoid singularities. On the one hand, both the Hartle- 
Hawking and Vilenkin wave functions can be taken gen- 
erally as non-zero at a = (being a minimum (Hartle- 
Hawking) or a Green's function (Vilenkin)), increasing 
or decreasing through the barrier. The singularity could 
be smoothed by means of some non-trivial measure in 
the inner product. On the other hand, quantum me- 
chanical situations with singularities occur and are dealt 
with, as with the solution of the Dirac equation for the 
ground state of the hydrogen atom which is singular when 
r — > 0. We have, however, herein chosen to follow the ap- 
proach presented in Ref. H^HMIHl to address the transi- 
tion from a dust phase towards an accelerating expansion 
stage by means of a quantum mechanical tunneling medi- 
ating the evolution from an initial Lorentzian (classically 
allowed) region to the current accelerating epoch. 

Second, it was also our objective to introduce a dis- 
cussion (and promote subsequent research) on the in- 
teresting similarities between the quantum cosmology of 
a specific brane-world models [3!j and FRW quantum 
cosmologies with a Chaplygin gas content. This was 
based on the feature that FRW models with a Chaplygin 
gas may constitute a phenomelogical tool to investigate 
brane-world models: both frameworks share a physical 



6 Let us mention that, in a different physical context, other quan- 
tization conditions (bearing formal similarities) were obtained in 
ref. l65l . when requiring normalization of the quantum states. 

7 The authors are grateful to C. Kiefer and A. Vilenkin for their 
comments on this issue. 



minisuperspace where two Lorentzian sectors are sepa- 
rated by a classically forbidden Euclidean region [55ll56| . 
This suggests a view that FRW models with a general- 
ized Chaplygin gas content may enable us to investigate 
specific aspects within brane- worlds. E.g., the quantum 
cosmological creation of the universe, by means of an 
Hamiltonian formulation, taking into account the pres- 
ence of potential symmetries and conserved charges that 
could be related with physical quantitites in brane-world 

physics [SE Eit Ea EE 13 Eli . 

This association between FRW cosmologies with a gen- 
eralized Chaplygin gas and some brane-world cosmolo- 
gies is strengthen when considering the next order of ap- 
proximation in Eq. (J5J and l|13|) (i.e., going beyond the 
Hamilton- Jacobi equation JTJJl): see Eq. (g7|) - lpT7|). This 
order of approximation was described through of a per- 
turbative expansion in the parameter A present in the 
Chaplygin gas energy density. Let us remind that a con- 
nection between the Chaplygin gas (a = 1) and brane- 
world models is that the brane tension is \[A [l£|- For 
the case of a generalized Chaplygin gas the "generalized" 
tension would be A 1 + a . A subsequent quantization con- 
dition involving A as well as a and B would be retrieved, 
depending on the boundary condition (Hartle-Hawking 
or Vilenkin) chosen. 

Let us now discuss the presence of the positive cosmo- 
logical constant, A, in our minisuperspace scenario. It 
may seem that its purpose is merely to induce the tun- 
nelling towards a DeSitter classical stage, as Eq. (fLUl sug- 
gests. Nevertheless, it is worthy to point the following. 
If we take A = in Eq. Jf)J) and depict V{a) for different 
values of a and with A <C 1, we get a physical situation 
as described in Fig. ^ (see Fig.[3|for a comparison), where 

V(a) 




a 



FIG. 3: This figure shows the behavior of the potential V(a), 
when the cosmological constant, A, is zero and the constant A 
is very small. The physical situation described in this case is 
similar to the one in Fig.Q i.e. two classically allowed regions 
disconnected and a Euclidean one. 

the behavior of an effective cosmological constant is in- 
duced through the parameter A. However, the analyt- 
ical quantum treatment of a minisuperspace with V(a) 
in Eq. ijfJJl with generic values of A, B and a is rather 
complex. The minisuperspace model described within 
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this paper with positive cosmological constant, A, is an- 
alytically feasible, providing a similar minisuperspace dy- 
namics through Eq. (fH5|l . The presence of A (pertinent 
for larger values of a) may be interpreted as the lowest 
order effect of the effective cosmological constant behav- 
ior induced in the generalized Chaplygin gas scenario. 
This provides explicit analytical expressions for the tran- 
sition amplitudes (from an initial Lorentzian (classically 
allowed) region to the current accelerating epoch through 
a tunneling process) to be computed. Improved expres- 
sions could be considered by taking next orders of ap- 
proximation, beyond the Hamilton- Jacobi Eq. (|10f) . as 
indicated in Sec. V. 

Another issue that is important to mention is the fol- 
lowing. In Eq. ijBJ we have a superposition of states, that 
may be considered from the point of view of decoherence 
(see Ref. [4fl | for an extensive and thorough review). In 
generic quantum cosmology, a full wave function could 

take the form W = S(„)C(„)e lS o" > x (n) where the are 
solutions to the Hamilton- Jacobi equation (treating the 
gravitational part) and the states x are matter states 

obeying i y 5*0™' V X (n) = # (n) X (n) , where is the 
corresponding matter Hamiltonian (see Ref. ^17]). Such 
expansion for '5 has some relevant non-classical features. 
On the one hand, it is a superposition of many WKB 
states. Only if those different branches are dynamically 
independent of each other we recover a background space- 
time without quantum mechanical interferences. On the 
other hand, restricting to a single component in ^ may 
not yet mean a classical space-time as the WKB function 
may spread over all configuration space. A decoherence 
process is needed thus prior to identify a classical behav- 



equation as in section V) in a broader perspective |5Cll52| . 
This would mean extracting the corresponding functional 
Schrodinger equation from the Wheeler-DeWitt equa- 
tion, together with quantum gravitational corrections to 
it (see Ref. |52| for more details). Moreover, the quantum 
cosmology of FRW models with a generalized Chaplygin 
gas could be investigated from the following perspective: 
taking into account the Chaplygin gas influence not as 
as additional term in the superpotential depending on a, 
but as a couple of independent conjugate dynamical vari- 
ables, which would imply a modification of the Wheeler- 
DeWitt equation. The guidelines of Ref. [28j could be of 
importance towards such Hamiltonian formulation. This 
will bring changes as the dimension of the configuration 
space will increase and the hyperbolic nature of the ki- 
netic term has to be considered. 



APPENDIX A: POSITIVE ROOTS IN CASE 

AB 2 < 4/9 

Let us write Vb(a) as 

Vo(a) = -Xaf(a) (Al) 

with 

f(a) = a 3 + b 2 a 2 + b 1 a 1 +b a°, (A2) 

where b 2 = 0, b\ = -A -1 , b = -f . Using Ref. [§]]], we 
define 



Concerning our case study, we focused our attention 
in the transition amplitudes for the universe to proceed 
from a classical initial region, through a quantum me- 
chanical tunneling sector, towards another classical re- 
gion where an accelerating stage takes place. We are tak- 
ing that different WKB components will decohere as dis- 
cussed in 0] an d considering solely a robust component 
(with respect to decoherence, in the sense introduced in 
Ref. |46l I47II49I ] ) of the original wave function. This com- 
ponent would induce a transition from a dust dominating 
stage towards an accelerated epoch as described in sec- 
tions I and II. A discussion on decoherence issues in our 
model (or other models where the Euclidean (classically 
forbidden) region contour is generically a closed curve) 
may be the subject of future research work. 

Finally, it may be of interest to indicate additional 
possible subsequent directions of research to explore. An 
interesting option is to inve stigate whether other brane- 
world models (see, e.g. Ref. [Hll00000llil) 
can admit a phcnomcnological description in terms of 
a generalized Chaplygin gas minisuperspace cosmology. 
This could be investigated by means of suitable canoni- 
cal transformations at the Hamiltonian formulation. An- 
other possibility is to pursue the analysis of the next 
orders of approximation (beyond the Hamilton- Jacobi 



A-\ r^-||^ g 3 + r2 = _ A - 3+ 9|J ; (A3) 



and then obtain for the case q 3 + r 2 < 0, i.e., AB 2 < |, 
the following roots for /(a) = 0: 



(A4) 
(A5) 
(A6) 



ai = 2 |a;| 3 cos ( - 



«2 



-2|x| 3 cos( - -- ) , 



03 = —2 |a;| 3 sin 



7T 9 

6^3 



In the last expression |a;| = A 2 and tan# = — ^3^ —. 
Moreover, we have that a\ — 0,3 > 0, together with 

ai > a 3 > > a 2 . (A7) 

Therefore, Vo(a) has three real roots determining two 
Lorentzian regions and one classically forbidden region 
for the range of a > 0. 
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